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Abstract. We present evidence that the mid infrared ( MIR , rest frame 5-30 /im ) is a good tracer of the 
total infrared luminosity, L(IR)(= L[8 — 1000 fim]), and star formation rate (SFR), of galaxies up to z ~ 1.3. 
We use deep MIR images from the Infrared Space Observatory (ISO) and the Spitzer Space Telescope in the 
Northern field of the Great Observatories Origins Deep Survey (GOODS-N) together with VLA radio data to 
compute three independant estimates of L(IR). The L(IR,MIR) derived from the observed 15 and/or 24 /im flux 
densities using a library of template SEDs, and L(IR, radio), derived from the radio (1.4 and/or 8.5 GHz) using 
the radio- far infrared correlation, agree with a 1-cr dispersion of 40%. We use the k-correction as a tool to probe 
different parts of the MIR spectral energy distribution (SED) of galaxies as a function of their redshift and find 
that on average distant galaxies present MIR SEDs very similar to local ones. However, in the redshift range 
z = 0.4-1.2, L(IR,24 /im) is in better agreement with L(IR, radio) than L(IR,15/im) by 20%, suggesting that the 
warm dust continuum is a better tracer of the SFR than the broad emission features due to polycyclic aromatic 
hydrocarbons (PAHs). We find marginal evidence for an evolution with redshift of the MIR SEDs: two thirds of 
the distant galaxies exhibit rest-frame MIR colors (L(12 /im)/L(7/im) and L(10 /im)/L(15 /im) luminosity ratios) 
below the median value measured for local galaxies. Possible explanations are examined but these results are not 
sufficient to constrain the physics of the emitting regions. If confirmed through direct spectroscopy and if it gets 
amplified at higher redshifts, such an effect should be considered when deriving cosmic star formation histories 
of dust-obscured galaxies. We compare three commonly used SED libraries which reproduce the color-luminosity 
correlations of local galaxies with our data and discuss possible refinements to the relative intensities of PAHs, 
warm dust continuum and silicate absorption. 
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1. Introduction 

Luminous infrared phases, during which galaxies radiate 
the bulk of their light in the infrared (IR) and are produc- 
ing stars at rates larger than 20 M yr _1 , play a dominant 
role in the cosmic star formation history (CSFH) and in 
the production of the cosmic infrared background (CIRB) 
JCharv k Elbadl200ll lElbaz et al.ll2002l Il999t IXu et alj 
l2003t iFranceschini et al.ll200ll 12003 lLagache et al.ll2004t 
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* Based on observations collected with the Spitzer Space 
Telescope (which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology under NASA contract 1407) 
and on observations with ISO (an ESA project with in- 
struments funded by ESA Member States (especially the PI 
countries: France, Germany, the Netherlands and the United 
Kingdom) with the participation of ISAS and NASA) 



ICharv et all l2004t IPanovich et all 12004ft . Those studies 
suggest that at least half of present-day stars were formed 
in "infrared luminous galaxies" , which we define here as 
the combination of luminous (LIRGs, 10 11 Lq < L(IR) < 
1O 12 L ) and ultra-luminous (ULIRGs, L{IR) > 10 12 Lq) 
infrared galaxies. Indeed, the comoving density of mid 
(15 & 24 /mi) and total IR luminosities (L(IR) = L[8 - 
1000 /im]) due to infrared luminous galaxies was about 70 
times larger at z ~ 1 t han it is today IjLeFloc'h et al.l2 005: 
ICharv fc Elba3l200lh . However, existing instruments do 
not allow direct measurement of L(IR) and one must 
rely on correlations between other wavelength ranges and 
L(IR) to derive its value as well as t he obscured star for- 
mation rate (SFR) that it quantifies l|Kennicuttl ll998) . 

L{IR) can be derived from a single measurement in 
the mid IR (MIR, rest frame 5-30 /im) for local galax- 
ies with a 68% uncertainty of 30-40% depending on 
the MIR wavelength (see Figs.l&2 from ICharv fc Elba3 



2 



Marcillac et al.: Mid infrared properties of distant infrared luminous galaxies 



(2001)). Similarly, the radio continuum exhibits an im- 
pressive correlation with L(IR) over four orders of mag- 
nitude w ith an uncertain ty of only 5% over nearly 2000 
galaxies (jYun et al .120011) . However, it must be noted that 
when limi ted to the 1 62 inf rared luminous galaxies of the 
sample of lYun et alJ l)200 lh . this dispersion rises to 46% 
which is comparable to the one observed between the mid- 
infrared and total IR luminosities. The origin of both cor- 
relations, linking the radio and MIR to L(IR), is not yet 
completely understood since several physical processes are 
generally advocated to explain the emission in both do- 
mains. 

It is generally acknowledged that the radio emission is 
produced by synchrotron radiation from relativistic elec- 
trons and free-free emission from HII regions. Most of the 
energy required for the acceleration of electrons is pro- 
duced by the supernova remnants of stars more massive 
than 8 M©, which also dominate the ionization in H II 
regions. However, whether this is enough to explain the 
extent of the correlation remains a matter of debate. 

The correlation between mid and far IR luminosities 
of galaxies is more direct since both are produced by the 
re-emission of UV photons by dust, but the physics of the 
emitting sources is quite complex. Above 60 /xm, the radia- 
tion of a galaxy is dominated by the grey body emission of 
"big dust grains" reaching an equilibrium temperature of 
20-50 K. At shorter wavelengths, the spectrum of a galaxy 
results from the combination of : 

- the continuum emission due to stochastically heated, 
warm dust grains, fluctuating around temperatures of a 
few hundred degrees. The carriers associated with this 
emissio n are commonly considered to be Very Small 
Grains foesert et al.lll990L VSG). 

- the PAH (polycyclic aromatic hydrocarbon) bands 
due to stochastically heated molecules producing broad 
features in emission at 3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 /itm. 

- the absorption features due to silicates which are cen- 
tered at 9.7 and 18 /im. 

Extragalactic surveys in the MIR, FIR, sub-millimeter 
and radio have all been used to trace star formation 
without being affected by dust extinction. However, with 
existing datasets and instrument sensitivities, the MIR 
is the most efficient wavelength to detect galaxies in 
the LIRG regime up to z ~ 2. Sub- millimeter surveys 
are limited to the mJy level where LIRGs are only de- 
tected up to z ~ 0.5 and with a bias towards objects 
with cooler dust temp eratures for a fixed IR luminosity 
1 Chapman et al . 2005). Deep radio surveys provided a pic- 
ture consistent wi th the one derived from the MIR but to 
a shallower depth jHaarsma et al.ll200ot lElbaz et al.ll2002l 
lElbaz fc Cesarskvl200.^ 7 As a result, it is crucial to deter- 
mine whether MIR properties of distant infrared luminous 
galaxies are similar to local ones and assess if the MIR can 
be used as a tracer of star formation in the distant uni- 
verse. This is the goal of the present paper. 

We will first compare the L(IR) derived indepen- 
dantly from three tracers, namely the radio continuum 
at 1.4 and/or 8.5 GHz, the 15/im flux density from 



ISOCAM dCesarskv et all I"l 996) onboard the Infrared 
Space Observatory ijKessler et al 1 ll996l ISO) and the 
24/im flux density fro m MIPS (Multib and Imaging 



Photometer for Spitzer iRieke e^ al.| |2Q04|) onboard the 
Spitzer Space Telescope ((Werner et all l2004() ■ The flux 
density and variation of the k-correction at the two MIR 
wavelengths are used to constrain the relative contri- 
butions of the three components of MIR light : VSGs, 
PAHs and Silicate abs orption. Using a sim ilar strategy, 
lElbaz et alJ l|2005h and iTeolitz et all l|2005h showed evi- 
dence that the MIR SED distant LIRGs did exhibit PAH 
features in emission. We will extend this work statistically 
by producing a direct comparison of the positions of local 
and distant galaxies in luminosity-luminosity diagrams at 
the same rest-frame wavelengths. This will allow us to test 
a possible variation of the relative roles played by the pre- 
vious three ingredients of MIR SEDs which would have 
strong implications on the derivation of the CSFH. 

While the derivation of L(IR) from a radio measure- 
ment depends only on the radio-far infrared (FIR) corre- 
lation and radio spectral index, which is well constrained 
for star forming galaxies, it is strongly dependant on the 
library of template SEDs that is chosen when using the 
MIR. Hence w e will compare three commo nly used S ED 
librar ies - CE llCharv fc Elbazll200lh. D H l(Dale fc HeloJ 
and LDP ijL^umcheetHini^l^ - which present 
some noticeable differences in terms of the relative con- 
tributions of the three major components of MIR light. 

We adopted a iJ =75 km s^ 1 Mpc" 1 , fl M = 0.3, fl A = 
0.7 cosmology throughout this paper. 



2. Sample 

The sample consists of 49 IR luminous galaxies de- 
tected with ISOCAM and MIPS, including a subsam- 
ple of 17 radio-detected galaxies, located in the Great 
Observatories Origins Deep Survey in the Northern hemi- 
sphere (GOODS-N) centered on the Hubble Deep Field 
North (HDFN). GOODS is a NASA Legacy Survey with 
Spitzer (Dickinson et al., in prep.) and with the Hubble 
Space Telescope (HST) Adv anced Camera for Surveys 
(ACS) llOiavalisco et al.ll2004l) . 

The size of the sample i s limited by the de pth of the 
ISOCAM catalog at 15 pm l)Aussel et al.lll999|) . This cat- 
alog is complete at the 90 % level down to 100 /iJy but a 
total of 95 sources are detected down to ^l OuJy among 
which 81% have a spe ctroscopic redshift l(Cohen et alJ 
| 2000HWirth et al.ll2004l) . We used a revised version of the 
lAussel et alJ l)1999l) catalog (H.Aussel, private communi- 
cation): 15/im sources were attributed an optical counter- 
parts from HST images and when two or more ISOCAM 
galaxies were closer than 12", the optical position of the 
HST sources was used to deblend the MIR sources. We 
rejected ISOCAM sources for which two potential optical 
counterpar ts were closer th an the FWHM of the ISOCAM 
PSF (4.6", IOkumurall998l) and kept only the sources with 
negligible residuals after deconvolution using the optical 
centers in the final catalog. After extracting those, the 
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catalog consists of 54 galaxies, all of which turned out 
to possess a MIPS counterpart. We checked all individual 
sources by eye in the MIPS-24 fim image and found that 
blending was not an issue for any of this sub-sample of 
galaxies. Finally, we removed the five active galactic nu- 
clei (AGN) that were identified from their X-ray emission 
^Fadda et all 120021: lAlexander et alJ 12002^ and the final 
sample consists of 49 galaxies. 

The radio data used in t his study c ame f rom VLA 
A-array observations done by I Richards! l)2000j) and new 
additional 26 hours of B-array data acquired by Morrison 
et al. (2006). The combined data sets and new reductions 
used here yielded an improved rms of 5.3uJy at the phase 
center of the reduced map and provided a significant frac- 
tion of new faint radio detections (Morrison et al. 2006) . 
The final sample is described in Tables 01 El 

The astrometry of the MIPS 24 ^m images is good to 
0.2" based on a cross correlation with the IRAC images. 
Fewer than 10 sources, none of which are relevant for this 
paper, are extended in the MIPS images which have a 
spatial resolution of 5.7" (PSF FWHM). Furthermore, al- 
most all 24 /im sources (down to ^20 /iJy, 5-er) are found 
to have an IRAC 3.6 /zm counterpart. 24 /im fluxes were 
determined by fitting a 24 /im PSF at the position of the 
IRAC sources, using a technique similar to DAOPHOT. 
The technique will be described in greater detail in the 
GOODS MIPS catalog paper (Chary et al. 2005, in prep.). 
This technique alleviates the effect of confusion noise and 
provides a direct association between the 24 /mi source 
and its shorter wavelength counterpart. 

The sample of local galaxies which will be compared to 
the distant luminous infrared galaxies consist of 154 galax- 
ies (z<0.1) with flux densities measured at 12, 25, 60 and 
100 /im with IR AS and at 7 and 15 M m wi th ISOCAM 
as described in ICharv fc Elbazl <|200ll) and lElbaz et alJ 
(2002). Only a subsample of these objects have 1.4 GHz 
radio flux densities. 

3. Total IR luminosity and SED libraries 

3.1. SED libraries and decomposition between VSGs 
and PAHs 

We briefly summarize in the following the origin of the 
three libraries of template SEDs - CE, DH, LDP - (see 
Fig. ^) and we refer the reader to the associated papers 
for more details. While the CE templates were built to 
reproduce the observed local luminosity-luminosity cor- 
relations from 7 to 850 /im, the DH ones were designed 
to fit the observed color-color relationships between local 
galaxies. The DH library consists of 64 templates sorted 
as a function of the fjjP/ft 00 ratio in the IRAS FIR 
bands. In order to compare this library to the CE and 
LDP ones, we attributed a total IR luminosity, L(IR) 1 
to each SED using Eq. ^ that we derived from the fit of 
the Ioq-[o\ ff Pl f,!°°l versus loq -\n\L(IR)] relationship (see 
Fig. 6a from ISoifer et all 17987) for objects in the Bright 
Galaxy Sample. The coefficients and dispersion that we 
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Fig. 1. Evolution of the SED templates with infrared 
luminosity in the 3-30 /im range top : Chary & Elbaz SED 
templates (2001) center : Dale & Helou SED templates 
(2002) bottom : Lagache et al (2004) 

used are given in Eq. ^ The resulting SEDs exhibit total 
IR luminosities in the range L(IR)= 2.1xl0 8 - 2.2xl0 14 
L Q (see Fig.P) and fit the local luminosity-luminosity cor- 
relations used by CE with a comparable quality. 

logic [fT/fT] = 0-128 x log 10 [L(IR)] - 1.611±°;^ (i) 

The shape of the 46 LDP starburst templates (L(IR)— 
3xl0 9 - 6xl0 12 L Q ) in the MIR present the smallest evo- 
lution with L(IR). They re sult from a revision by LDP 
of the iLaeache et alJ l)2003|) templates where the contin- 
uum in the 12-30 /im has been modified to allow a better 
fit of the MIPS-24 /im number counts. Anothe r 46 n or- 
mal templates were also built bv ILaeache et alJ l|2003|) to 
reproduce the population of spiral galaxies which radiate 
more than half of their energy in the optical. Since in- 
frared luminous galaxies do not fall in this category, we 
did not use these SEDs for the estimation of L(LIR,MIR). 
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L,r (L ) 

Fig. 3. Evolution of the Lpah /(Lpah+Lvsg) with L(IR) at 7/im (up) and 12 /urn (down). Sample of galaxies : 0: 
SBS 0335 (0.03 Z ) 1: IIZw 40 (0.2 Z Q ), 2: NGC 5253 (0.2 Z Q ), 3 : NGC 1140 (0.4 Z ), 4 : NGC 520, 5 : M 51, 
6 : NGC 3256, 7 : Cen A, 8 : M 82, 9 : IC 342 (central starforming regio n), 10 : NGC 6946 11 : NGC 253, 12 : 
Circinus, 13 : NGC 1365, 14 : NGC 6240, 15 : M 8 3 16 : NGC 4038/3 9 llMirabel et allll998l overlapping regio n 
of the Antennae galaxie s), 17 : nucleus of NGC 40 38 l|Mirabel et alJll998(l . 18 : Arp 220 ()Charmandaris et al.lll999j) . 
19 : IRAS 23128-5919 l|Charmandaris et al.ll200^) . Galaxies number 4 to 19 have metallicities between 0.6 and 2 Zq 
while galaxies number to 3 have metallicity < 0.4 Zq as indicated within parenthesis. 



However, we will present their location in the MIR color- 
color diagrams of Sect. 03 

We have used the me thod of IVerstraete et aL | ll200lh 
and iGalliano et all l|2005f) to decompose each SED from 
all three libraries into the sum of broad emission lines at 
5.3, 5.7, 6.2, 7.7, 8.6, 11.3, 12.0, 12.7, 13.6 and 14.3mn due 
to PAHs and warm dust continuum due to VSGs, here- 
after called the PAH and VSG contributions. Lorentzian 
profiles w ere adopted for the emission features (as pro - 
posed by iBoulanger et all ll99St IVerstraete et all 1200 lj) . 
with line widths and intensities as free parameters. The 
VSG continuum was modeled with a modified black body, 
oc B V {T) x v° , (3 — 1 being the emissivity; the temper- 
ature, T, and the intensity of this component are free to 
vary, in order to minimise the \ 2 of the fit. The decom- 
position of the DH SED with L{IR) = 9xlO n L Q is pre- 
sented in Fig. |21 This technique provides excellent fits for 
all SEDs regardless of luminosity and without introducing 



a silicate feature in absorption at 9.7 /mi. The presence of a 
trough around 10 /im that is visible mostly in the CE SEDs 
below 5xlO n L Q (Fig. is perfectly reproduced by the 
depression between the right wing of the 8.6 pm, feature 
and the left wing of the 11.3 [im one. It has been shown 
by various authors that the MIR SED of local star form- 
ing galaxies can be reproduced in the absence of silicate 
absorption at 9.7 fim which is difficult to s eparate from 
PAH bands when the latter are intense (e.g. Sturm et al. 
|2000|) . PAHs are always strong in the three SED libraries 
both around 7 /im and 12 /xm, except for the most lumi- 
nous galaxies in the CE library, but this is also where no 
trough is noticeable (see Fig. 01 . 

This decomposition allows to follow the relative con- 
tribution of each component of the SEDs, i.e. PAHs and 
VSGs, as a function of wavelength for all SEDs and to 
study its evolution with L(IR). We computed this evo- 
lution with L(IR) in two passbands centered on the two 
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Fig. 2. Demonstration of the decomposition method on 
the template galaxy with L(IR) = 9 x 10 11 L®, in the DH 
library. The spectrum (thin line) is fitted with the linear 
combination of a hot grain continuum (bold line), PAH 
bands (dashed lines) and ionic lines (dotted lines). 



ties (se e iGalliano et all (|2003t 120051) : lEngelbracht et alJ 
l|2005a|) . 



3.2. Derivation of the total infrared luminosity 
(8-1000 urn) 

The correlations between the FIR (40-120 /im, £fir) or 
total IR (8-1000 /im, L{IR)) luminosity with the radio 
and MIR can be used to derive independent estimates of 
L(IR). We use L(IR, radio), L(IR,15) and L{IR,2A) for 
the infrared luminosity derived from the radio, ISOCAM- 
15 and MIPS-24 passbands respectively. Lfir and L{IR) 
are relat ed by Eq . El comp uted from the Bright Galaxy 
Sa mple llSoifer et alJ | l98^) where we use the definitions 
of iHelou et alJ l|l988j) for L F ir and ISanders k. Mirabel 
l|l99rt for TATR). 

L IR = 1.91(±0.17) x L(FIR). (2) 

For the radio-FIR cor relation, we used a q-paramete r 
of 2.34 as defined in Eo.l3l l|Condonll992tlYu"n et al.l2001^ : 



strongest PAH features at 7 and 12 /im. These two pass- 
bands correspond to the rest-frame central wavelenths 
probed by the ISOCAM-15/xm and MIPS-24 /im filters at 
z = 1. The two components for each passband are here- 
after called L£ AH , Lj SG and i^ H , LYi G . The DH and 
LDP libraries present a similar behavior (see Fig. |3J), i.e. 
very little variation of the contribution of PAHs to these 
bands over the whole range of total IR luminosities: PAHs 
are responsible for ^90 % (~50 %) of the light in the 7 /im 
(12 /im) band. 

On the contrary, while the contribution of PAHs re- 
mains larger than 80% in the 7 /im band for the CE li- 
brary, it becomes smaller than that of the VSGs in the 12 
/im band L(IR) — 5xlO n L Q , and reaches 20% for the 
highest luminosities. A direct comparison with a sample 
of local galaxies for which both ISOCAM spectroscopy 
and IRAS broadband fluxes are available suggests that 
the behavior of the DH and LDP templates is more repre- 
sentative of the 7-12 /xm color for local galaxies although 
we do see a trend for the contribution of PAHs to de- 
crease in the ULIRG regime (see galaxies 18 and 19 in 
Fig. . These local galaxies are comprised of local sp irals 
ijRoussel et al.ll200]T) . starbursts jLaurent et alJl2000l) and 
dwarf galaxies l|Madden et aLlEoORbf) observed with the 
Circular Variable Filter of ISOCAM. Note that the LDP 
templates provide a ratio of the 12 over 7 /xm flux densities 
which does not vary with luminosity, contrary to what is 
observed for local galaxies, and is twice smaller than the 
observed average ratio (see Sect.EJ). 

In Fig.0 galaxies number 4 to 19 have metallicities be- 
tween 0.6 and 2 Z®, typical of relatively massive galaxies 
at < z < 1. Four galaxies (numbers to 3 in Fig.|3J), have 
metallicity lower than 0.4 Z®. Their properties may be re- 
flective of primeval galaxies, for which local SED libraries 
might not be valid in the derivation of total IR luminosi- 



_ , ( L(FIR)(W) 1 \ 

q l0gw [S.75 x 10" (Hz) L^GKziWHz- 1 )) ■ { ) 

The rest- frame 1.4 GHz monochromatic luminosity was 
derived from the observed 1.4 or 8.5 GHz flux densities 
assuming a radio s pectral slope for star forming galaxies 
of a = 0.8 ± 0.15 , ilYun et al.ll200ll S v cx v~ a ). 

To derive L(IR) from the MI R (15 and /or 24 /xm) , 
we used the strategy advocated in ICharv fc Elba! ll200lh . 
Each one of the template SEDs is redshifted to the dis- 
tance of the observed source. Then, the SED whose MIR 
flux density is the closest to the observed one is used to 
derive L(IR) after it is normalized to the exact observed 
flux density. 

4. Global consistency of the 15/iim, 24/im and 
radio star formation indicators 

The 15 and 24 /im passbands probe different sources of 
MIR emission at the median redshift (z ~ 0.85) of the 
sample considered here. At this distance, the 15 /im band 
samples the PAH emission centered around 7.7 /im while 
the 24 /*m band measures the combined emission from 
VSGs and the 11.3 and 12.7 /xm bands. However, both 
bandpasses provide consistent values for L(IR) (see Fig.0} 
with a median ratio of 1 and an rms of 40 % when using 
the DH library (LIRGs/ULIRGs only, i.e. third line of 
Table 1). This dispersion is the same as the MIR versus 
L(IR) one for local galaxies which suggests that the MIR 
SED of galaxies did not strongly evolve with redshift up 
to z ~ 1. 

The 18 galaxies of the sample for which both MIR 
and radio data are available are located at redshifts be- 
tween 0.3 < z < 1.3. The total IR luminosity derived from 
24 /xm and the radio are consistent within 40 % (rms de- 
fined as the 68% confidence level, see Fig. Et, i n agree- 
ment with previous studies l|Elbaz et alJ 120021 iGarrettl 
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Fig. 4. a : LIR\§ as a function of LIR24 for the SED of CE; b : Evolution of the ratio LIR15/LIR24 with respect to 
LIR24 for the SED of DH; c : LIR U as a function of LIR 2 a for the SED of CE; d : Evolution of the ratio LIR 15 /LIR 24: 
with respect to LIR24 for the SED of DH; the size of the circles is proportionnal to z. 









L{IR, 15)/L(/.R,24) 
CE DH LDP 


all galaxies 
L(IR,24)< 10 11 
L(IR,24)> 10 11 


Lq 
L 


0-841^ 
0A7 + _° o f 2 
0.85±°;« 


1 9+0.8 

1 - A -0.7 

i-3l l :g 


1-Oto-l 


0.5l° ;l 



Table 1. Median and 68% dispersion (around the median) of the ratio L(IR, 15)/L(IR, 24) estimated with the SED 
of CE, DH and LDP. 



l2002t iGruppioni etaD 120031: lAppleton et all l2004|) . The 

median ratio of L(IR,MIR) over L(IR, radio) is ~ 1 (see 
Table |2J) but the dispersion around this value depends 
on the MIR wavelength that is used and on the li- 
brary of template SEDs. Note that the 40 % dispersion 
of L(IR,MIR)/L(iR,radio) is comparable to the 46% dis- 
persion for L(FIR)/ L(IR n radio) found for t he 162 local IR 
lumin ous galaxies of lYun et al.l l(200l|) and lReddv fc Yur] 
(2004). Hence, the agreement between MIR and radio de- 
rived total IR luminosities at z ~ 1 is at least as good 
as in the local universe, which suggests that the SED of 
galaxies in the IR were not strongly different at Z ~ 1 
from what it is today. 



The 24 /im passband gives the lowest dispersion, with 
an rms of ~ 30 and 40% for the CE and DH libraries 
respectively. At 15 /im, the dispersion becomes as large as 
100 % for the CE library and 60 % for DH. The very good 
agreement of L(IR,24) with L(IR, radio) supports the idea 
that the radio remains a robust SFR indicator up to z ~ 
1.3. This in turns implies that the 15 /im passband is a 
less robust tracer of L(IR) than 24 /jm. Since the median 
redshift of those galaxies is z ~ 0.85, this suggests that 
the warm dust continuum is a better tracer of L(IR) than 
PAHs. 

A comparison between the results obtained with the 
CE and DH libraries shows that systematic effects arc ob- 
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Fig. 5. Comparison of the total IR luminosity (8-1000 /im) derived from the radio (from the 1.4 and/or 8.5 GHz 
bands) and the MIR at 15 /im (open circles) and 24 /im (filled circles). Figs a,c and Figs b,d were produced using 
the CE and DH libraries of template SEDs respectively, to derive I/t^mir- The solid line in Figs. a, b is the one-to- 
correlation between the two L(IR). The solid and dashed lines in Figs c,d represent the median (thick lines) and 1-a 
(68 %) dispersion (thin line) around this median for the 24 (solid) and 15 /jm (dashed) data respectively. 

z L(IR, MIR) /L(IR, radio) 

15 /im 24 /im 

CE DH LDP CE DH LDP 



10 11 <L(IR) < 10 12 L Q 0.85 t°oi 1-2 ±o.t 1.0 t° 6 s 0.6±°; 3 3 1.0 ±°; 3 2 1.1 ±° A 3 l.3f 



Table 2. Median and 68% dispersion around the median of the ratio L(IR, MIR) / L(IR, radio) estimated with the 
template SEDs of CE, DH and LDP. 



served for L(IR,24) > 10 11 L Q . We find that the L(IR,15) 
derived from the CE library is 30% larger than L(IR,24) 
and that the ratio of both luminosities presents a disper- 
sion between 60 and 90% (see Table 1). The major dif- 
ference between the CE and DH libraries is the progres- 
sive dilution of PAHs by the VSG continuum as a func- 
tion of increasing L(IR) for the CE SEDs (see Fig. Pi If 
this dilution were to be excessive or if PAHs were to be 
stronger than predicted by this library then one would ex- 
pect the CE SEDs to overestimate the L(IR) derived from 



a flux density measured around the peak of PAHs emis- 
sion around 8 /im. This would explain the excess ratio of 
L(IR,15) over L(IR,24). And indeed, when compared to 
an independant tracer of L(IR), the radio, we find that 
L(IR,15)/L(IR, radio) is 20% larger than one, with a dis- 
persion larger than with the DH library (see Table 0) . 
When computed using the LDP library, L(IR,24) is twice 
larger than L(IR,15) (see Table 1). The origin of this dis- 
crepancy is opposite to the one with CE, i.e. the strength 
of PAHs around A ~ 8 /im (in the rest-frame of the z ~ 
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0.85 galaxies) is too strong for these luminosities. Indeed 
the L(IR,MIR) predicted with the LDP library exhibit a 
systematic offset with respect to L(IR,radio) (see TablcEJ). 
To avoid confusion, we did not plot those points in Fig. [5] 
This discrepancy results from the 7.7 fim PAH feature 
being too strong and the continuum plus PAH emission 
around 12 /xm being too faint in their SED library. Note 
that LDP designed their library in order to better repro- 
duce the 24 /jm galaxy counts observed with Spitzer. 

Finally, we note that the three galaxies with redshifts 
close to z ~ 0.5 - for which L(IR, radio) < 3xlO n L Q - 
exhibit a L{I R, 15) / L(IR, radio) ratio lower than 1 both 
for CE and DH (Fig.|5J. At z ~ 0.5, the 15 /im passband 
probes the 10 /im trough in the rest-frame which might be 
underestimated by those templates. The number of galax- 
ies at this redshift for which both MIR and radio data are 
available is too small in the present sample to derive any 
strong conclusion. However, we will address this point by 
comparing the MIR sample to local galaxies in the next 
section. 

5. Testing the redshift dependence of 
mid-infrared SEDs 

In order to check whether distant and local luminous IR 
galaxies exhibit consistent MIR properties, we have com- 
pared the evolution of the observed 24 /zm over 15 /im lu- 
minosity ratio with redshift to prototypical starbursts and 
active galactic nuclei (fig. EH) and to SEDs such as CE and 
DH (fig Et) ■ We have also selected two redshift bins (z ~ 
0.6 and 1) such that the 15 and 24 /im passbands probe 
rest-frame wavelengths for which data on local galaxies 
are already available. Local and distant galaxies are rep- 
resented with filled triangles and circles respectively in the 
Figs. 03 El We distinguish three domains in luminosity for 
"normal" galaxies, "LIRGs" and "ULIRGS", correspond- 
ing to [10 10 -10 n ], [ 10 n -10 12 ] and greater than 10 12 L 
range. In both figures, the median ratio of the luminosi- 
ties for local and distant galaxies are marked with dark 
solid triangles and circles respectively. The lcr error bars 
were computed by the quadratic sum of a poissonian error 
and an error computed to include 68 % of the galaxies in 
each bin. We will center the discussion on the normal and 
LIRG populations due to a lack of objects outside of these 
domains. We will start by presenting the redshift evolu- 
tion of the observed 15 over 24 /im color as a test for the 
presence of PAHs in distant galaxies. 

5.1. Testing the presence of PAHs in distant galaxies 
using the 15 over 24\im ratio 

In Fig. EH, we compare the location of the GOODS-N 
galaxies (circles) in the 24 over 15 /im luminosity ratio ver- 
sus redshift space to several pro totypical local SEDs. This 
figure reproduces the Fig. 3 of lElbaz et alJ l)2005|) from 
which the square points are taken. The closest starburst 
M82 (L(IR) = 4 x 10 10 L Q ), the merging system of the 
Antennae (NGC 4038-39; LIRG, L{IR) = 1.4 x 10 11 L Q , 




0.0 I i i i I i i i I i I i I i i i I i I i I i i i I i 

0.0 0.2 0.4 0.6 0.8 1.0 1,2 

redshift 

Fig. 6. vL x ^ m /vL^ m as a function of the redshift. 
Black and grey circles : our distant sampl e of galaxies. 
Black and grey squares : data presented in lElbaz et alJ 
l|2005h . Red circles : : median i/L^ m /vL^™ in each 
redshift bin with lcr error bars. 



IMirabel et afl lll998hL the highly obscured LIRGs NGC 
4418 dSpoon et all200ll). the classic, b ut atypical, ULIRG, 
Arp 220 l|Charrri;mdaris et al.1 1999) and the Compton 
thick Sevfert II NGC 1068 l|Le Floc'h et all200l[) are com- 
pared to our distant LIRGs. 

The templates can be divided into three families: 

1. Starbursts. M 82 and th e overlapping regio n 
of the Ante n nae galaxies IMirabel et alJ 
IVieroux et all l|l996h : the SEDs of M82 and of the 
overlapping region of the Antennae galaxies (Mirabel 
et al. 1998, Vigroux et al. 1996) globally reproduce 
the behavior of the observed data but tend to re- 
spectively over- and under-estimate the L15/L24 ra- 
tio above z=0.5. This could be explained either by the 
warm dust continuum being too faint or the silicate 
absorption feature being too weak. 

2. Heavily obscured starbursts: Arp 220, NGC 4418 and 
the nucleus of NGC 4038. 
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NGC 4418 is either a strongly embedded AGN or a 
compact starburst (Evans et al, 2003). The MIR spec- 
trum shows strong silicate absorption without strong 
evidence for PAH features ( Spoon et al, 2001). 
The MIR spectrum of Arp 220 shows characteristics of 
a highly obscured continuum with silicate absorption 
(like NGC4418) and a starburst component character- 
ized by PAH emission (like M82) (Spoon et al, 2004). 
All three cases predict a low value for L15 / L24 at z<0.9 
and a high value at z >1.1. This is probably due to a 
combination of strong silicate absorption and contin- 
uum emission from warm dust. 
3. Active galactic nucleus (AGN). NGC 1068 harbors a 
central black hole with a circumnuclear starburst. 75 
% of the MIR emission is dominated by the AGN emis- 
sion (Le Floc'h et al, 2001). The MIR spectrum is 
dominated by a strong continuum, without PAH. Its 
L15/L24 ratio is too weak for 68 % of the sample due 
to its flat SED and the absence of PAH features in 
emission. 

In Fig. EJd, we have indicated three bins of red- 
shift (0—0.5, 0.5—0.9, 0.9—1.3) which correspond approx- 
imately to three bins in luminosity because of selection 
effects: [3 x 10 8 - 3 x 10 11 ], [3 x 10 9 - 7 x 10 11 ] and 
[10 10 — 3 x 10 12 ] Lq. In each bin, the red dot is associated 
to the median value of the data with an error bar including 
68 % of the sample. The position of the red dots confirms 
the conclusion of Elbaz et al. (2005) that the L15/L24 ra- 
tio increases from z ~ 0.2 to 0.7 due to the entrance of 
the 7/im PAH emission feature into the ISOCAM-15/im 
band. 

The template SEDs from CE and DH are represented 
for three luminosities, L(IR) = 10 10 8 (thin line) , 10 11 ' 4 
(thick line) and 10 124 (very thick line) L . The global 
behaviour of the data is well represented by both families 
of synthetic templates. However, for the CE templates, 
the z = 0.5—0.9 bin is better represented by SEDs less 
luminous than LIRGs while the observed galaxies are pre- 
dominantly LIRGs. This confirms the previous statement 
that in the CE templates, the PAHs emission in IR lumi- 
nous galaxies is underestimated by about 15-20%, while 
this is not the case of the DH templates. 

5.2. L\2 versus Lj 

At z ~ 1, the 15 and 24 /im observations measure rest- 
frame wavelengths of approximately 7.5 and 12 /im, re- 
spectively. We will designate the rest-frame luminosities 
at these wavelengths as L7 and L\i. Lj measures the 
"PAH bump" centered on the 7.7 /zm dominant PAH fea- 
ture, while L12 measures a combination of the 11.3 and 
12.7/mi PAH bands with the warm dust continuum due 
to VSGs, as well as some contribution from the silicate fea- 
ture in absorption, if any. We may compare L7 and Lyi for 
the high-redshift galaxies to data for local (z ~ 0) galax- 
ies observed with ISOCAM-LW2, centered at 6.75 /im, 
and IRAS 12 fim respectively. Fig. shows the rest- 




10 8 10 9 10'° 10" 

L, 2 (L S ) 

Fig. 8. Comparison between the rest-frame 7 and 12 urn 
luminosities of local (filled triangles) and z ~ 1 ( 0.8 
< z < 1.25 ,big filled circles : fl5 > 150 fiJy, little 
filled circles : fl5 < 150 fiJy) galaxies. L7 and L12 
correspond to vL v at these central wavelengths for the 
local- ISOCAM-6.75 /im and IRAS-12 /jm passbands - and 
z rsj 1 - rest-frame bands corresponding to the observed 
ISOCAM-6.75 fim and MIPS-24/im - galaxies. Two do- 
mains of luminosities are defined by vertical lines (see text) 
where the local and distant galaxies are compared. (a)Lj 
as a function of L12. The plain line is a polynomial of de- 
gree 2 fit to the local sample alone, (b) L12/L7 ratio as 
a function of L\i- The median ratio for local and distant 
galaxies are indicated with dark triangles and circles re- 
spectively, in each luminosity bin. The lines indicate the 
prediction for the CE (plain line), DH (dash dot line), 
LDP cold (thick dashed line) and starburst (dashed 
line) templates. 



frame wavelengths observed by the 15 and 24 /jm bands 
at z ~ 1 as grey shaded regions, and compares them to 
the ISOCAM-LW2 and IRAS 12/xm bandpasses for lo- 
cal galaxies. The comparison between local and distant 
galaxies at these wavelengths is presented in Fig. El Since 
the redshifted ISOCAM and MIPS filters do not exactly 
match the observed ISOCAM and IRAS filters used for 
local galaxies, we computed a color correction for each in- 
dividual galaxy, taking into account its proper redshift, 
using the SED libraries. This color correction factor is al- 
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Fig. 7. Relative positions of the ISOCAM-15 /xm and MIPS-24 /xm bands for galaxies located at z ~ 1 (a) and 0.6 (b) 
with respect to the observed ISOCAM-LW2, LW3 and IRAS-12 /xm bands, represented in the plots is /„ in arbitrary 
units as a function of A. Origin of the template SEDs at 3 different luminosities (L(IR) — 10 10 L Q) "normal", 10 11 
L Q , "LIRG" and 10 12 L , "ULIRG"): CE (Solid line), DH (Dash-dotted line), LDP (Dashed line: Starburst 
template, thick dashed line "cold" template). 



ways smaller than 20 % and the choice of the CE or DH 
library makes less than a 5 % difference. 

We extracted a subsample of 116 galaxies from the 
local sample, all of which have flux measurements at both 
7 and 12 /xm. The distant sample was made of 26 distant 
galaxies with a flux at 15 and 24 /xm and located in the 0.8- 
1.25 redshift range in the redshift range 0.8 < z < 1.25, 
with median redshift z = 0.90. 

Distant galaxies follow the same trend as local galaxies 
in their L7 versus L12 correlation (Fig. Eh.)- This shows 
that at least in the MIR domain, the SED of galaxies 
have not evolved strongly between z ~ and z ~ 1. 
We do find an indication of marginal evolution: nearly 
all z ~ 1 galaxies lie above the polynomial fit of degree 
2 to the local sample (solid line in Fig. |H^). We com- 
puted in Fig. [Ht> the median L12/L7 ratio for local (filled 
black triangle) and distant (filled black circle) LIRGs. 
Distant galaxies exhibit a lower L12/-L7 ratio than local 
ones in both luminosity bins (LIRGs- normal galaxies), 
with an apparently stronger difference in the lower lumi- 
nosity bin. However, this bin consists of only six galaxies, 
four of which are fainter than 75 /xJy and therefore lie be- 
low the 60 % comp leteness regime of the ISOCAM image 
llAussel et al Il999l) . The risk that the galaxies detected at 
this flux density level might not be representative of the 
average galaxy properties (Malmquist bias in the sense 
of favoring objects with strong 15 /xm flux densities and 
low L12/L7) is therefore not negligible. Hence we consider 
that our sample is not robust against biases outside of the 
LIRG regime. 

In the LIRG luminosity bin, all but two of the galax- 
ies fall above the 90 % completeness limit of 100 /xJy. The 
median L12/L7 ratio for local galaxies lies above the 1-cr 
dispersion (including error bars) around the median ratio 
for distant galaxies, i.e. two thirds of the distant galaxies 



present a LnjL-j ratio lower than the median value for 
local galaxies. Although, this effect needs to be consoli- 
dated by a larger sample, we believe that it is important 
to consider the possibility that we are already seeing some 
indication of evolution that will become stronger at larger 
redshifts. Possible explanations for this evolution are dis- 
cussed in Sect. El 

The L12/L7 versus -L12 trend is globally fitted by the 
CE and DH libraries but not the LDP one, which sys- 
tematically underestimates L12/L7 over all luminosities. 
This observation confirms the excess PAH emission, al- 
ready discussed in Sect.QJ for this library. Finally, we note 
that the CE library presents a rapid increase of L^/Lj 
with increasing luminosity in the IR luminous regime and 
that the L(IR) derived using this library might be overes- 
timated by up to 50% for distant ULIRGs. This can be 
explained by the same reasons discussed in the previous 
section, i.e. the weak emission from PAHs around 8 /xm in 
this library. 

5.3. Liq versus L15 

At z ~ 0.6, the 15 and 24 /xm observations measure rest- 
frame wavelengths of approximately 9.4 and 15 /xm, re- 
spectively. We will designate these rest-frame luminosities 
as L10 and L15. While both Liq and L15 contain emission 
from the VSG warm dust continuum, only Liq is affected 
by the silicate absorption feature. Hence, the flux ratio 
of these two bands helps to quantify the strength of the 
silicate absorption. We may compare Lxs f° r high-redshift 
galaxies to data for local (z ~ 0) galaxies observed with 
ISOCAM-LW3, centered at 15/xm. There is no direct ob- 
servation of Liq for local galaxies, but we may derive it 
from a combination of observations with ISOCAM-LW2 
and IRAS-12/iin as described below. We will designate 



Marcillac et al.: Mid infrared properties of distant infrared luminous galaxies 



11 




\ 2 - 



LisCU) 



Fig. 9. Comparison of the rest-frame 10 and 15 /im lu- 
minosities of local (filled triangles) and z ~ 0.6 (0.4 < z 
< 0.8, big filled circles : fl5 > 150 fjjy, little filled 
circles : fl5 < 150 /iJy) galaxies. L15 corresponds to 
vL v at this central wavelength for the ISOCAM-15 /im 
and MIPS-24 fan passbands for local and z ~ 0.6 galaxies 
respectively. L\§ corresponds to the rest-frame emission 
measured in the ISOCAM-15 /im band at z ~ 0.6 and 
was derived from a combination of the IRAS-12/im and 
ISOCAM-15/im bands for local galaxies (see Eq. 0] (a) 
L10 as a function of L15. The plain line is a polynomial of 
degree 2 fit to the local sample alone (here equivalent to a 
linear fit), (b) L10/L15 as a function of L15. The median 
ratio for local and distant galaxies are respectively indi- 
cated with dark triangles and circles in both luminosity 
bins. The lines indicate the prediction for the CE (plain 
line), DH (dash dot line), LDP cold (thick dashed 
line) and starburst (dashed line) templates. 



this synthetic bandpass as "LW3-0.6". We illustrate the 
local and redshifted bandpasses for the z ~ 0.6 sample in 
Fig. Cb- The comparison between local and distant galax- 
ies at these wavelengths is presented in Fig. El 



We have compared 115 local galaxies that have mea- 
sured 12 and 15 /xm flux densities to 35 distant galaxies 
at z ~ 0.6 (0.4 < z < 0.8). In order to probe the wave- 
length range associated with the 10/um feature, we con- 
structed a synthetic passband centered at 10 /mi using a 
combination of the ISOCAM-LW3 (12-18 /im) and IRAS- 
12 (im (8-15 /im) filters. The passband of the synthetic 
filter was chosen to reproduce the rest-frame wavelength 
range probed by the ISOCAM-LW3 filter for a galaxy lo- 
cated at z ~ 0.6, i.e. centered at 9.4 /im and therefore 
close to the 10 /im feature . We found that a combina- 
tion of the ISOCAM-LW3 and IRAS- 12 /mi filters could 
provide an estimate of the synthetic passband (hereafter 
called LW3-0.6) flux density with an accuracy of 20 %. We 
optimized this combination for CE and DH templates in- 
dependantly but found that the difference between both 
computations is only 10% (see Eq.0]computed for galax- 
ies which L(IR) > 4.7 x 10 10 L Q ). The LW3-0.6 filter is 
represented with a line in fig. At z = 0.6, the MIPS- 
24 /im passband is roughly equivalent to the rest-frame 
ISOCAM-LW3 filter; however, we calculated a filter cor- 
rection to take into account the k-correction and the dif- 
ference of filter shapes. The correction is estimated at 10 % 
when using the DH templates and 30 % for CE. 



T CE 

n LW3-0.6 



1.64 
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0.71 



-0.70 
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Here again, as in Sect. 15.21 distant galaxies (light grey 
circles) follow the trend traced by local galaxies (filled tri- 
angles) , hence confirming that no strong evidence is found 
for an evolution of the MIR SED of galaxies with redshift 
(see Fig. EJ). Again, distant galaxies exhibit a systemati- 
cally lower Lio/£i5 ratio with two thirds of the distant 
LIRGs showing a lower ratio than local LIRGs. This can- 
not be explained by a selection effect since in this lumi- 
nosity range, all galaxies exhibit ISOCAM flux densities 
above the 90 % completeness limit. Possible explanations 
are discussed in Sect. 

6. Discussion and conclusions 

We have presented a comparison of the mid infrared (MIR) 
and radio properties of local and distant galaxies (up to 
z ~ 1.3) with particular emphasis on IR luminous galax- 
ies, i.e. LIRGs and ULIRGs. Our goal was to assess the 
possibility of deriving a total IR luminosity, L(IR), corre- 
sponding to a star formation rate (SFR), on the sole basis 
of a measured flux density in the MIR for those galax- 
ies. We derived three L(IR) estimates using independantly 
the measured ISOCAM-15 fim, MIPS-24 fan and VLA-1.4 
GHz flux densities and found that all three luminosities 
were consistent at the 40 % confidence level, i.e. an uncer- 
tainty similar to the one measured for the local MIR-FIR 
and radio-FIR correlations for IR luminous galaxies. This 
first result, obtained for a sample of 49 galaxies with red- 
shifts up to z ~ 1.3, and a sub-sample of 18 galaxies de- 
tected in the radio, suggests that the MIR can be used as 
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a tracer of the SFR of distant actively star forming galax- 
ies. However, the L(IR) computed from 15 (im shows a 
10-20 % larger scatter than the one computed from 24 /im 
in comparison to the radio derived L(IR). This result sug- 
gests that the warm dust continuum and neutral PAH is a 
better tracer of the SFR than is the ionized PAH emission 
at rest-frame 7.7 /im. 

The redshift evolution of the 24 over 15 /zm colors show 
a bump around z ~ 0.7 as one would expect from the 
entrance of PAH emission centered around 7.7 /im into the 
15 um. passband. This beh avior reinforces prev ious studies 
bv lElbaz et afl l)2005|) and lTeplitz et all (|2005|) which were 
also interpreted as evidence for the presence of PAHs in 
distant IR luminous galaxies. 

We attempted to probe the same rest-frame wave- 
length range in local galaxies and high redshift galax- 
ies. We selected galaxies within two redshift bins cen- 
tered at z ~ 0.6 and 1, where the observed 15 /im (24 /im) 
passband corresponds to the rest-frame 10 and 7/im (15 
and 12 /im) respectively. In both redshift bins, the dis- 
tant galaxies fall in the same region as local galaxies of 
equivalent luminosities and provides an extension of lo- 
cal correlations to higher luminosities. We find marginal 
evidence of evolution with redshift of the MIR SEDs: 
two thirds of the distant galaxies exhibit rest-frame MIR 
colors (L(12/im)/L(7/im) and L(10/im)/L(15/zm) lumi- 
nosity ratios) below the median value measured for local 
galaxies. Possible explanations are examined below. It is 
indeed to be expected that an evolution should take place 
if not at the z ~ 1 level, then at higher redshifts. The first 
cause that one may consider is the role played by metallic- 
ity. It was recently showed bv lEneelbracht etah (2005b) 
and lMadden et all l|2005a|) that the PAH over VSG con- 
tinuum ratio was increasing as a function of metallicity, 
with a trend for galaxies with metallicities below log(Z)= 
12-Mo<?(0/H)=8.2 to exhibit nearly no PAH contribution. 
iLiang et alJ l)2004() found the gas metallicity of LIRGs 
ranging between half-solar to solar (log(Z)=8.63-8.93) 
at similar redshifts. Hence, the influence of metallicity is 
not expected to be strong in our sample. Furthermore, it 
would be responsible for a stronger L12/L7 for the dis- 
tant galaxies instead of the lower values that we tend to 
find. We measure a L^/L? ratio for the distant LIRGs 
which is similar to the one found for "normal" nearby 
galaxies. This may suggest that distant LIRGs, which are 
probably more gas rich than local ones, could also be trig- 
gered by less violent causes, i.e. minor mergers, gas in- 
fall or even spiral waves instead of the major mergers 
foun d in local galaxi es, as suggested by their morphologies 
fsee lBell et al.ll2005lk We note that the dust temperature 
measured in the far infrared also de pends on the geome- 
try of the regions of star formation ijChanial et al J 12005: 
IWang fc Heloulll992l) . The lower L10/L15 ratio may result 
from a deeper 9.7 /im silicate absorption. Another possibil- 
ity would be that the level of ionization of PAHs in distant 
galaxies is higher: more ionised PAHs tend to pr esent a 
lower L12/L7 fe.g. lDartois fc D'Hendecourtlll997l . Other 
causes might also play a role such as a different grain size 



distribution but it is not possible to disentangle between 
all those possibilities with existing data. Extended sam- 
ples of objects at higher redshifts with Spitzer 16 /im and 
24 /im imaging as well as mid- infrared spectroscopy will 
soon be available which might help discriminate between 
the various effects. 

We teste d three commonly used librari es of template 
SEDs- CE (iGharv fc Elbad'hoPl!). DH llDale fc HeloiJ 



120021) and LDP l|Lagache et al.ll20ol - in which the rel- 
ative strengths of broad emission lines from PAHs with 
respect to warm dust continuum due to VSGs increases 
from CE to DH and to LDP. 

The LDP library predicts values for L(IR) which are 
40 % too low from the measured 15 /im flux density and 
30 % too large from the 24 /tm one, with respect to the ra- 
dio ones. They also exhibit the largest dispersion and in- 
consistency between L(IR,15) and L(IR,24). This is mostly 
due to enhanced 8 /tm and decreased the continuum be- 
tween 11-32 /zm and and to a lesser extent, due to a lack 
of evolution of MIR SEDs as a function of L(IR). 

The CE and DH libraries both provide robust fits to 
both local and distant galaxies but we find evidence for 
more PAH emission in IR luminous galaxies than pre- 
dicted in the CE SEDs. As a result, the CE library would 
overestimate L(IR) for distant galaxies by 20-30%. If, as 
suggested above, distant galaxies do exhibit stronger PAH 
emission than local ones, then both libraries need to be re- 
vised and the CE library would overestimate L(IR) by up 
to 50% for ULIRGs at z ~ 1. Since LIRGs dominate the 
star formation rate density per unit comoving volume at 
least up to z ~ 1.5, this evolution of MIR SEDs does not 
significantly affect the derivation of the cosmic star forma- 
tion history (CSFH). We have ch ecked this by replacing 
the CE SEDs by the DH ones in the lCharv fc Elbazl l|200lK 
model and found two CSFH which are consistent within 
30%. 

The true accuracy with which FIR luminosities can 
be determined at high redshift are being tested for 
the most luminous objects with Spitzer, and will be 
pushed to fainter luminosities with the Herschel satel- 
lite scheduled for launch 2007, and the Atacama Large 
Millimeter Array, ALMA. Direct MIR spectroscopy with 
the InfraRed Spectrograph onboard Spitzer will also pro- 
vide very useful constraints if pushed to detect distant 
LIRGs. 
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